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Abstract 

We consider an special dynamics of a quantum walk (QW) on a line. Initially, 
the walker localized at the origin of the line with arbitrary chirality, evolves to an 
asymptotic stationary state. In this stationary state a measurement is performed and 
the state resulting from this measurement is used to start a second QW evolution to 
achieve a second asymptotic stationary state. In previous works, we developed the 
thermodynamics associated with the entanglement between the coin and position 
degrees of freedom in the QW. Here we study the application of the first and second 
laws of thermodynamics to the process between the two stationary states mentioned 
above. We show that: i) the entropy change has upper and lower bounds that are 
obtained analytically as a function of the initial conditions, ii) the energy change is 
associated to a heat-transfer process. 
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1 Introduction 


Quantum walks (QWs) constitute the quantum analogue of classical random 
walks [lj and also the quantum version of cellular automata |2j. They have 
been intensively investigated, especially in connection with quantum informa¬ 
tion science . As in the classical case, QWs have been proposed 

as elements to design quantum algorithms [TUilllll2llT3] and more recently 
it has been shown that they can be used as a universal model for quantum 
computation [1JJT5]. 


We have been investigating [1611171118] the asymptotic behavior of the QW on a 
line, focusing on the probability distribution of chirality independently of po¬ 
sition. We showed that this distribution has a stationary long-time limit that 
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depends on the initial conditions and that it is possible to define a thermo¬ 


dynamic equilibrium between the degrees of freedom of position and chirality 



[T2I12 0112 11122] . For this equilibrium state we have introduced a temperature 

concept for an unitary closed system. 

On the other hand, the fundamental lower bound of the thermodynamic energy 
cost of information processing has been a topic of active research [231124] . On 
average, the minimum amount of work required to erase 1 bit of information 
from a memory is kT In 2 [25]. In the last decades developments in nano-science 
have enabled the direct measurement of such minuscule amounts of work for 
small non-equilibrium thermodynamic systems [26]. At the same time recent 
advances in technology have opened the possibility of building useful quantum 
computing devices [27] . Therefore, it seems essential to identify bounds on 
the thermodynamic energy cost of information processing [28] for these new 
quantum devices. 

In the present paper we study the relationship between the QW thermody¬ 
namics and information processing. In particular, we show that it is possible 
to apply the thermodynamic laws to the QW dynamics after a measurement 
process. We obtain the upper and lower bound for the asymptotic change 
of the entanglement entropy. Our result may be thought as complementary 
to the results presented in Refs. j2 M] where the information content and 
thermodynamic variables are treated on an equal footing. 

The paper is organized as follows. In the next section the usual QW on a line 
is presented, in the third section the system dynamics with measurement is 
developed, in the fourth section the entropy change between the asymptotic 
stationary states is studied, in the fifth section the laws of thermodynamics 
are applied to the same process. Finally, in the last section we draw some 
conclusions. 


2 QW on a line 

The composite Hilbert space of the QW is the tensor product T-Lt <8> H± where 
T-Lt is the Hilbert space associated to the motion on a line and T~L± is the 
chirality (or coin) Hilbert space. In this composite space the walker moves, at 
discrete time steps t 6N, along a one-dimensional lattice of sites k G Z. The 
direction of motion depends on the chirality states, either right or left. The 
wave vector can be expressed as the spinor 



( 1 ) 
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where the upper (lower) component is associated to the left (right) chirality. 
The QW is ruled by a unitary map whose standard form is 


a k (t + 1) = a k+1 (t) cos 6 + b k+1 (t) sin 9, 
b k {t + 1) = a k -\(t ) sin# - b k _i{t) cos 0, 


( 2 ) 


where 6 G [0,7t/2] is a parameter dehning the bias of the coin toss, ffere we 
take 0 = j for an unbiased or Hadamard coin. The probability of finding the 
walker at (fc, t ) is 

P(k,t) = \a k (t)\ 2 + \b k (t)\ 2 . (3) 

The global left and right chirality probabilities are defined as 


p L(t)= Y M*)l 2 , 

k=—oo 


Pr(*)= E \h(t)\ 2 , 

k=—oo 

with Pn(t) + Pl( t) = 1 and the interference term is dehned as 

OO 

Q(t) = Y. a k (t)b* k (t). 


( 4 ) 


( 5 ) 


k=— oo 


In the generic case Q{t ) together with P L (f) and -Pr(£) are time depend func¬ 
tions that have long-time limiting values [IB] which are determined both by 
the initial conditions and by the map in Eq.([2]). The relation between the ini¬ 
tial condition and the asymptotic distributions has also been recently explored 
in Ref. [33]. Let us call the mentioned limits as 


hr = 


Pr = 


Qo = 


lim P L (t) 
t —> oo 
lim P R (t) 
t —> oo 
lim Q{t) 
t — y oo 


= ii + iv, 


( 6 ) 


where n and v are respectively the real and imaginary part of Qo. The following 
relations are verified [16] between LR, II R and Q 0 


n L ,- 

n h = 1 


/i, 




(7) 
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It is important to emphasize that the asymptotic behavior in Eq.([6|) is deter¬ 
mined by the interference term Q 0 that only depends on the initial conditions. 


The initial condition the walker localized at the origin with arbitrary chirality 
will play a central roll in our analytic treatment. Then 


*i(0)> 


cos(7/2) 

^expip sin( 7 / 2 ) 


| 0 >, 


( 8 ) 


where 7 G [ 0 , 7 r] and 99 G [0, 27t] define a point on the unit Bloch sphere. We 
obtain [IB] 


Q 0 


-(1- y=) COS 7 + sin 7 (cos^j + i\/2sinp) 

2 v 2 L 


( 9 ) 


and also II L and II # using Eq.©. 


3 Dynamical evolution and measurement 

3.1 First step 


We consider first the QW evolution starting from the initial condition given 
by Eq. (l 8 ]l and determine the asymptotic density matrix. The quantum density 
matrix is defined as 

p(f) = |T(f))<T(f)|, ( 10 ) 

substituting Eq.([T]) into Eq. ffTOjl we have 


P(t) = J2 

k,k' 


f a k (t)a* k ,(t) a k (t)b* k ,(t) 
K a* k (t)b k f(t) b k (t)b* k ,(t) 


I k)(k' 


( 11 ) 


where a k {t) and b k {t) depend also on the initial conditions 7 and 99 . 

Note that when t —> 00 the limits of a k (t) and b k {t) are not defined because, in 
general, they have an oscillatory asymptotic behavior [ 3 ]; however the limits 
given by Eq.(j6]) are always well defined. I 11 the following we call a k and b*. to 
the values of a k (t) and b k (t) respectively, evaluated at times large enough so 
that the asymptotic limit, Eq. flB]) . is essentially attained. 


Let us define the asymptotic reduced density matrix as p\ c = limtr(p) = 


4 






lim 5Dz(Z|/?|/) for t —> oo. This matrix takes the following shape 


Pic 


V Qo 

yQo 


( 12 ) 


3.2 Second step 


When the asymptotic density matrix is attained a measurement of position 
and chirality is performed. Then the wave function collapses into one element 
of the set of eigenvectors of the measurement operator (|fc)|±)}, where 


l+> = 



(13) 


After the measurement, the density matrix has the form 

OO OO 

Pi— E \*k\ 2 \k+)(k + \+ ^ \hk\ 2 \k—)(k — \, (14) 

k =—oo k =—oo 


where |/c±) = |fc)|±). Eq. flTTl) can be written in matrix form as 

/ 


pi = E 


V 


|a fc | 2 0 
0 Ibfcl 


\k)(k\. 


(15) 


Equation ([IT]) sets a new initial condition and we will study the asymptotic 
evolution of the new density matrix. Essentially we must study the evolution 
of the densities |/c±)(fc±|. We know, from the QW dynamics, that if the initial 
condition is given by Eq. ([8|) then the density matrix is given by Eq. ([TT]) . We 
can use these equations as a recipe to obtain the unknown evolutions, that is, 
when we choose a particular initial condition given by | k±) we have 


U{t)\k±){k±\U\t) = 


E 

n,n' 


' a n ±(t, k)a* n ,±(t, k ) a n ±(t, k)b* n ,±{t, k) 
K a n *±(t, k)b n >±(t, k) b n ±(t , k)b* n ,±(t , k) 


n)(n '\, 


(16) 


where U(t) is the QW evolution operator and a n ±(t, k) and b n ±(t , k) verify the 
map Eq.(|2]), that is a n ±(t,k ) and b n ±(t,k ) are equivalent to some functions 
a n (t) and b n (t). The expressions a n ±(t, k) and b n ±(t,k) show explicitly their 
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dependence both on the chirality, (|+), |—)}, as with, k, the walker’s initial 
position on a line. Therefore, the new quantum density matrix is 


Pit) 


U{t) Pl U\t) 


E 

n,n' 


* Aiiin, n') 

> 

to 

s' 

s. 


\n){n 

yA 21 {n,n') 

A 22 (n,n') J 


where 


A u (n,n') = [|a fc | 2 a n+ (t, k)a* n , + (t, k) 

k 

+ |bft| 2 a n -(t,k)a’ n ,_(t,k )], 


( 17 ) 


(18) 


A 22 (n,n') [\ak\ 2 b n+ (t, k)b* n , + (t, k) 

k 

+ |b fc | 2 b n _(t, k)b* n ,_(t,k)] , 


(19) 


A 12 (n,n') = Y^ [|a fc | 2 a n+ (t, k)b* n , + (t, k ) 

k 

+ |b(:| 2 a„_(i, k)b‘ n ,_(t, *,')] . 


( 20 ) 


A 21 (n,n') =A‘ I2 (n,n'). 


( 21 ) 


We again point out that the probability density has an asymptotic limit for 
long times. We now call the values of a n ±(t, k) and b n ±(t, k) evaluated at such 
times a n ±{k) and b n± (£;) respectively. We want to calculate in this limit the 
reduced density matrix, namely p 2c = limtr(p) = lim t —> oo. Using 
the density matrix given by Eq. flTTlh we have 


Pic 


lim J2 


OO 

l=—oo 


An(l, l ) 2l 12 (/, l) \ 

^Al 2 (l, l ) A 22 (l, l) J 


t —)■ oo 


( 22 ) 


where 


li m E£Loo^n(M) 
t —> oo 


= El a fcl 2 El a «+( fc )l 2 

k l 


+ El b ^l 2 El a *-( fc )l 2 > 

k l 


(23) 
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(24) 


limEKLoo^MU) 

t —y OO 


=J2\ a k\ 2 J2 a i+( k ) h h( k ) 

k l 


+J2\ h k\ 2 J2 a i-( k ) h *-( k )’ 

k l 


limEg-oo/W) 

t —> OO 


= zl a fcl 2 El b «+( fc )l 2 

fc z 


+ El b fcl 2 El b *-( fc )l 2 - 

fc z 


According to Eqs.dU), (|5]) and ([6]), we define 

n L± = Sla, ± (*)| 2 , 

Z 


(25) 


(26) 


n fl± ^Elb.±W| 2 , (27) 

Z 

Oo± = E a <±w b ;±w- < 28 ) 

z 

To obtain the explicit shape of n L± , and Qo±; we can use again Eqs.([7l 
[9]) where the initial condition is, the walker localized at the position k with 
chirality |±) respectively, i.e. 


S'2±(0)) = |±)|fc). 


(29) 


Note that the principal difference between Eq. (JH]) and Eq. ()29l) is in the walker’s 
initial position, that is, in Eq. fl29jl k is arbitrary and in Eq.Q k = 0. However, 
even with this difference, Eq.Q continues to be valid for the calculation of Qo± 
because the original map, Eq. (J2]) , is invariant under translations and therefore 
this magnitude is independent of k. For the same reason Eqs. (l26|) and (1271) are 
independent of k , and then Eq. ()22l) reduces to 

P2c = hi lPl + TIrPr, (30) 


where 


Pl 


n l+ Qo+ 
yQo+ n H+ 


(31) 


Pr 


^n L _ q 0 _ 
^Qo- n«_ 


(32) 
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and U L and are given by Eqs. (0 [9]). Moreover, using the initial condition 
of Eq. (l29|) in Eq.((9]) and Eq.([7]) it is straightforward to obtain 


Qo+ 2 2V2' 


(33) 

Note that pl and pr, Eqs.d2HE2]), are the asymptotic densities that correspond 
to initial conditions associated with the eigenvalues |+) and |—) respectively. 
Therefore, Eq. (|30l) has a natural interpretation, after the measurement the 
asymptotic reduced density matrix is the weighted average of the two densities 
associated with the two possible values of the chirality. 


n,L+ — n fi _ — i — 


n fi+ — ru_ — 


2 V 2 


2 V 2 ' 


4 Entropy change 


The unitary evolution of the QW generates entanglement between the coin and 
position degrees of freedom which can be quantihed through the associated 
von Neumann entropy. This entropy of entanglement is defined by the reduced 
density operator 

S(Pc) = -« tr(p c log p c ), (34) 

where k denotes a proportionality constant, the Boltzmann constant. 

For a pure state the minimum entropy S(p c ) = 0 is attained. The maximum 
entropy (or minimum purity) is to be found for the broadest possible proba¬ 
bility distribution, the equipartition over all pure states. 

Equation (j5Tj) can be expressed as a function of the eigenvalues of p c , A + and 
A_ 

_ = — A + logA + — A_ log A_. (35) 

K 

Using Eqs. ffl2ll . (j30l) . (13B and ()32]i we can calculate the entanglement en¬ 
tropies for the four stationary densities introduced in the previous section. 
The eigenvalues of the density operator p ic are 

A 1± = \ ± V 2p 2 + (36) 

A 2± = i±H(V2-l). 


and those of p 2c are 


(37) 







The operators pl and p R have the same eigenvalues and they are 


11/3 

A LR± ”2 T 2Y~ 71' 

Therefore, the four entanglement entropies are expressed as 
S(pic) 

sCptc) 

K 

S(PL ) 

hi 


(38) 


K 


- A,, log Ai + — Ai_ log Ai_, 

(39) 

- A 2+ log A 2+ - A 2 _ log A 2 _, 

(40) 

— A lr+ log A L r+ — A lr ^ log A LR _, 

(41) 

S(Pi) « 0.139 . 

(42) 


hi 


We now show that S(p 2c ) has a non obvious upper bound, using a known 
theorem that plays an important role in many applications of quantum in¬ 
formation theory [34] . Applying the mentioned theorem to the entropy of the 
mixture of quantum states given by Eq. (j30j) . the following upper bound is 
obtained 

*5*2 = S(p L ) -n(U L logU L + U R logU R ) > S(p 2c ). (43) 

The entropies S(pi c ), S(p 2c ) and the expression S 2 only depend on the inter¬ 
ference term Qq. 

Note that from Eq. (J9]) we have |Qo| < |(1 — < 1, then we approximate 

the log functions in Eqs. (l39]h (140|) and (l43lh using the first few terms of their 
Taylor series to calculate the entropy change between the two asymptotic 
stationary states, 


^ 2u 2 + 2(2^2 - l)p 2 . (44) 

n 

The distance between S(p 2c ) and its upper bound can be approximated by 

& - K gW ^ 2 . 

K K 

Combining Eqs. (J43J), (144)1 and (145|) . it is easy to obtain the following upper 
bound for S(p 2c ) — S(pi c ) 

Ji = S(p R ) + 2 n(p 2 + v 2 ) > S(p 2c ) — S(pi c ). (46) 


When Q 0 vanishes, Eq. (l44]i shows that S(pi c ) = S(p 2c ) and these entropies 
take their maximum value. Moreover in this case ( Q 0 = 0) the QW dynamics 
can be described as a classical Markovian process it has a Markovian 
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behavior both before and after the measurement. The initial conditions for 
this behavior are 7 = — 7t/4, 7t/ 4 and = 0, 7 r respectively. When Q 0 7 ^ 0, 
the measurement process essentially determines an entropy increase for the 
system, S(p 2c ) - S(p lc ) > 0 . 



Fig. 1. Dimensionless entropy change as a function of |Qo | 2 = P 2 f° r the case tp = 0. 
In thick line the entropy change given by Eq. (1441) . In dashed line the entropy change 
given by Eq. (l45l) . 

The entropy change for a real Qo is depicted in Fig. [I] The calculation was 
made using the entropy definition Eqs. (13911401 [4TL 1421) , therefore the curves 
show that the lineal approximation proposed in Eqs. (14411451) is excellent. 


5 Thermodynamic laws. 


I 11 previous works [191201211122] . we studied the behavior of the chirality dis¬ 
tribution and we introduced the temperature concept for an unitary closed 
system. In this theory one considers the system associated with the chirality 
degrees of freedom and described by p c , interacting through its entanglement 
with the position degrees of freedom, the lattice, as equivalent to the thermal 
contact with a bath. Therefore, in equilibrium 

[H c , Pel = 0, (47) 

should be satisfied, where H c is the interaction Hamiltonian between the chi¬ 
rality and the lattice. In the QW case, the explicit shape of H c is unknown for 
us, however we know that the eigenvalues of H c are independent of the initial 
condition of the wave function, they only depend on the unitary evolution. In 
contrast, the eigenvalues A± depend on the initial conditions and the corre¬ 
sponding eigenvalues of H c . We call {|<f>±)} the set of eigenfunctions common 
to the density matrix and the Hamiltonian, then in this basis the operators p c 
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and H c are both diagonal. Moreover, since only the relative difference between 
energy eigenvalues has physical significance, we denote this set of eigenvalues 
by ±e; they may be interpreted as the possible values of the entanglement en¬ 
ergy. This interpretation agrees with the fact that A± is the probability that 
the system is in the eigenstate |$±). 

The precise dependence between A± and ±e is determined by the type of 
ensemble we construct. We propose that our equilibrium state corresponds to 
a quantum canonical ensemble. To this end we set 

A± = (48) 

where Z is the partition function of the system, that is 

Z = e~P £ + e fle , (49) 

and the parameter /3 can be put into correspondence with an entanglement 
temperature 

k(3 Klog(A + /A_) ^ 

The entanglement temperature, Eq. ()50lh can take any finite, infinite, positive 
or negative value. Note that when the energy of a system is bounded from 
above there is no compelling reason to exclude the possibility of negative 
temperatures. 


In this statistical mechanic frame it is possible to define the internal energy 
of entanglement between the coin and position degrees of freedom 

U(p) = eA + — eA_. (51) 

The variations of temperature, AT and internal energy, AU between the two 
asymptotic states can be calculatedas as functions of the interference term Q 0 . 
Using Eqs. (l36j) and ()37li together with Eqs. ()50]) and (j5Tj) we have 

_ _1_1_ C52) 

e log (Ai+/Ai_) log (A 2 +/A 2 _) 

1 1 

i/2p 2 + v 2 (^/2 — l'j \p\ 

^ = (V2 - l) \p\ - (53) 

Equation f)52|l shows that AT = 0 is only possible if A 1+ A 2 _ = A 1 _A 2+ , 
and this implies that Qo = 0, then the two asymptotic states have the same 
temperature and in this sense we can to think of an “isothermal process”. As 
the temperature concept makes sense only in an equilibrium state it is clear 
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that between these two asymptotic states the temperature is not defined. 
Additionally Eg. (|53|) implies that AU < 0 and vanishes when Q 0 = 0. 

The first law of thermodynamics is now applied to the evolution between the 
two asymptotic stationary states expressed as 

A U=Q + W, (54) 

where Q is the QW heat absorbed during the measurement process and W is 
the work made over the QW. Thermodynamic work is defined to be measurable 
solely from the knowledge of external constraint variables. In this system the 
only parametric dependence of the thermodynamical functions is with the 
temperature, because the energy levels {—e,e} are maintained constant, then 
W = 0. Therefore the first law is reduced to 

AU = Q. (55) 

From Eq. (l55l) we conclude that: i) The change of internal energy is due to 
the heat delivered during the measurement process, ii) The heat behaves as 
a state function, then the real process can be substituted by a quasi-static 
process between the two asymptotic stationary states (characterized by their 
temperatures), where it is possible to defined the infinitesimal dQ. 

AS i 
0.1 
0.01 
0.001 
0.0001 
IE-005 

IE-006 
IE-007 
IE-008 
IE-009 
IE-010 


Fig. 2. (Color on line) The dimensionless entropy change (in log scale) as a function 
of the dimensionless interference term Q o = [x + iv. The medium (green) surface is 
the change of entropy given by Eo. (l44l) . The upper and lower (red) surfaces are the 
upper and lower bounds Ji and J 2 given by Eq. fHBj) and Eq. (f59j) respectively. 

In general, for an irreversible process, the second law of thermodynamics is 
expressed as 

*»/%■ 
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(56) 




The question may be posed if this law is satisfied for the entropy, temperature 
and heat that were defined in our system. The answer to this question must 
be positive because the starting point of our theory has been to postulate 
the canonical distribution for the reduced density matrix therefore the ther¬ 
modynamic laws are obeyed. We do not know the temperature dependence 
with the absorbed heat in order to compute the integral in Eq. (l56|) . However 
we can obtain a bound for this integral in order to verify the second law for 
the entanglement between the coin and position degrees of freedom. We have 
|Tj < \T 2 \, where T x and T 2 are the temperatures of the stationary states 
before and after the measurement. Then 

■h = -|r < / ^ < Sfec) - Siftu). (57) 


Therefore, using Eq. (j55l) we can propose the second law for the entanglement 
entropy 

S(p2c) - S{pi c ) > J 2 = > 0. (58) 

J- 2 

The lower bound J 2 can be expressed as a function of the interference term 


Qo, Eqs. ([36U37U50U531 [551 [58D 
h 


K 


(V2-l)M - (V2 - l) PI 


(59) 


In Fig. [21 we show the entropy change together with its bounds. The entropy 
calculations were made using the definitions Eqs. (1391 [40]) . We conclude that 
the second law of thermodynamic in the form of Eq. f[58|) is satisfied by the 
theory developed in the present paper and additionally we show the correctness 
of the upper bound given by equation Eq. (l46l) . 


6 Conclusion 


In previous works we developed the thermodynamics associated with the en¬ 
tanglement between the coin and position degrees of freedom. Here we consider 
a special dynamics of a QW on a line. Initially, the walker localized at the 
origin with arbitrary chirality, evolves to an asymptotic stationary state, then 
a measurement is performed and the state resulting from this measurement is 
the initial condition for a second QW dynamics to achieve a second asymptotic 
stationary state. 

We have studied the first and second laws of thermodynamics in the process 
between the two stationary states mentioned before. These asymptotic sta¬ 
tionary states only depends on the initial conditions through the interference 
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term of the initial wave function. We show that the change of entropy has up¬ 
per and lower bounds and they are obtained analytically as a function of the 
initial conditions. We have also shown that the measurement process changes 
the energy and this change is associated to a heat transference process. 

Moreover, we prove that, if the interference term vanishes the thermodynamics 
functions of the asymptotic stationary states do not change. 


I acknowledge stimulating discussions with Victor Micenmacher, and the sup¬ 
port from PEDECIBA (Programa de Desarrollo de las Ciencias Basicas) and 
ANII (Agencia Nacional de Investigation e Innovation) (FCE-2-211-1-6281, 
Uruguay), and the CAPES-UdelaR collaboration program. 
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